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The correlated density appears in many physical systems ranging from dense interacting 



gases up to Fermi liquids which develop a coherent state at low temperatures, the su- 
perconductivity. One consequence of the correlated density is the Bernoulli potential in 
superconductors which compensates forces from dielectric currents. This Bernoulli po- 
tential allows to access material parameters. Though within the surface potential these 
, contributions are largely canceled, the bulk measurements with NMR can access this 

■ potential. Recent experiments are explained and new ones suggested. The underlying 

qq ' quantum statistical theory in nonequilibrium is the nonlocal kinetic theory developed 

' earlier. 

o 

1 
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1. Introduction 



The theory of dense quantum gases distinguishes the number of free particles and 
the number of particles in a correlated state, e.g. in bound states. Therefore any 
kinetic theory should result in a balance equation which in equilibrium leads to 
q \ the Beth-Uhlenbeck equation of state 1-3 where the density of free and correlated 

particles are distinguished n = n,f + n COIT (nf). In contrast to that the Landau 
theory of quasiparticles describes dense interacting Fermi systems only in terms of 
the quasiparticle density. This is supported by the Luttinger theorem 4 which states 
■ that the Fermi momentum for any correlated Fermi system in the ground state is 

exclusively determined by the total number of particles and there is no shift due 
to correlations. This apparent contradiction between the two approaches is solved 
within the nonlocal kinetic theory which unifies both pictures and provides them 
as limiting cases of a more general treatment. 5 How is this possible? First of all an 
ideal Fermi liquid obeys the Luttinger theorem and any correlated density has to 
vanish in the ground state. This is in a sense a tautology. 6 But one has to observe 
that most Fermi systems do not turn to ideal Fermi liquids at low temperatures. 
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Electrons in metals and nucleons, e.g. develop a coherent state, the superconducting 
state. There is a small fraction of correlated density in the superconducting state. 
Such correlated density shifts the Fermi momentum and consequently the chemical 
potential /i. Since the electrochemical potential has to remain constant there must 
be an electrostatic potential (j> compensating this shift according to fx + eip = 0. 
This electrostatic potential turns out to be of the form of a Bernoulli potential in 
superconductors. If this potential could be measured we have a proof of the concept 
of correlated density. 

This overview discusses the theoretical treatments of the Bernoulli potential 
and the experimental realizations to measure it. First we describe shortly how the 
correlated density appears in superconducting states and how it is related to the 
electrostatic potential. Then we review the history of this Bernoulli potential and 
explain why surface measurements have not been successful in the past. We found 
the reason in the Budd-Vannimcnus theorem which explains why large-scale cancel- 
lations appear on the surface of superconductors. Then we calculate what potential 
can be measured deeper in the bulk by calculating the density and magnetic profile 
of superconductors. This allows us to explain a recent NMR experiment and to pre- 
dict inhomogencous fields to be measured above vortices in type-II superconductors. 
Finally we conclude with a short description of the nonlocal kinetic theory which 
provides the description of the correlated density in nonequilibrium. 

2. The concept of correlated density 

2.1. Correlated density within BCS approach 

In superconductors the Wigner distribution function has a two-part structure 7 

fw = \ (l + |) f(E) + f(-E) = \-^ tanh \ E (1) 

where £ = e p — /i is the free-particle energy e p minus the chemical potential /i. The 
quasiparticle energy E = ^/£ 2 + A 2 describes the influence of the superconducting 
gap A on the excitation spectrum of the superconducting state. The Fermi-Dirac 
distribution is f(x) = l/(e l3x + 1) with the inverse temperature (3 = l/ksT. There 
are two branches of energies. For states inside the Fermi sphere, £ < 0, the negative 
branch given by the second term dominates. For the states above the Fermi sphere 
the role of the contributions is reversed. Close to the Fermi sphere both contributions 
are comparable. 

The density n is obtained by the momentum integral over (1) and introducing 
the density of states h(£) = 2 J2 p 27rr5(£ — e p ), one has 

oo 

»=/i^+oQ-A tanh ^)- (2) 

-A* 

Here we account for a possible electrostatic potential ip and the velocity v of su- 
perconducting electrons by p, — fi — eip — mv 2 /2. For a vanishing gap we obtain 



February 6, 2008 22:30 



WSPC - Proceedings Trim Size: 9.75in x 6.5in corrdense 



the corresponding density n n of normal electrons with the chemical potential p by 
n n = n(A = 0). The difference 

^corr — ^ (3) 

describes the correlated density. Let us discuss it first in the ground state. There 
the normal density turns into 

n„ = 2^e(Al-e p )«no-(e^+y« 2 )^ (4) 
p 

where we have expanded p in first order around the Fermi energy and no describes 
the number of particles with no motion and no electrostatic potential. The correlated 
density (3) splits into two parts in the zero-temperature limit of (2) 

o 



2 { 27r y/Z 2 + A 2 2 J 2tt vp w ^2 + A 2 V 7 

which vanishes for vanishing gap. Since the gap is only nonzero in the vicinity of 
the Fermi level given by the Debye frequency ujd we can restrict the integration to 
the ±w_D-range. Expanding the density of states for £ < too we obtain finally 



_ dh A 2 



ln ^ + ^ + 1 



9/iA 2 (2cj d 



for 3> A in the last step. 

Since the total system should stay neutral we expect n — n and the two con- 
tributions, n n — no according to (4) and n CO rr of (6), should cancel. Therefore the 
required electrostatic potential must read 

m , d In h A 2 , / 2lo d \ 

e<P = -^-v 2 + — — In \ . (7) 

^2 dp 2 V\/eA/ V 7 

We see that it has the form of a Bernoulli potential which compensates the contri- 
bution of the supercurrent velocity and the associated inertial and Lorentz forces 
and has a part directly linked to a material parameter namely the gap. The latter 
is called thermodynamic corrections. This internally developed Bernoulli potential 
and field keep the charge neutrality and compensate the inhomogeneities due to 
diamagnetic currents. 

The great hope is now to measure the Bernoulli potential in order to access 
directly the gap parameter. 



2.2. Bernoulli -potential 

A way to discuss the Bernoulli potential closely connected to experimental facts is 
the London theory. The London condition 



mv = — eA 



(8) 
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fixes the trajectories of superconducting electrons due to the transverse vector po- 
tential. This condition explains the Meissner effect and contains both London equa- 
tions. One can consider it as the simplest way to account for basic electrodynamic 
features of superconductors. The time derivative of this London condition gives the 
force 

dA ( 1 \ 

rav = -e— - e(vV) A = e(E + vxB) + Vle<^+ -mv 2 j (9) 

exerted on the electrons. Comparing this with Newton's equation of motion mv = 
e(E + v x B) + F s where the unknown force on the superconducting pair is denoted 
with F s one obtains 

\7&p = F s - X7^mv 2 . (10) 

The different expressions in literature differ by the different treatment of the un- 
known force F s . London neglected this force consistent with the London condition. 
This corresponds to the picture of hydrodynamics of a charged ideal gas and it re- 
sults in eip ~ — ^mv 2 . The diamagnetic current is maintained by the electrostatic po- 
tential. 8 Sorkin 9 suggested that the free energy responsible for superconductors must 
contribute, too, which yields eip ~ -|mw 2 — ^f-. As a further step van Vijfeiken 
and Staas 10 considered the force resulting from the interaction between electrons 
and condensate acting on electrons to keep them at rest F„ + eE = 0. Taking into 
account the balance of forces n n F n +n s F s = as well as that eE = — e\7ip the force 
becomes F s = — ^eVip. Introducing this into (10) one obtains that the Bernoulli 
potential is reduced by the fraction of superconducting electrons eip = —^■^rnv 2 
called quasiparticle screening or fountain effect. Rickayzen 11 finally started from 
the kinetic energy density fkin = n s ^mv 2 of the condensate which determines the 
change in the chemical potential 

d dn s 1 2 n s 1 2 n n d\nT c l 2 

eip = -u = -—fkin = — ■K-^ mv = n mv + 4 -mv A . (11) 

an on 2 n 2 n a Inn 2 

This expression provides, besides the quasiparticle screening, also thermodynamic 
corrections. This offers the experimental possibility to access directly material pa- 
rameters similar to the BCS expression (7). 

The experimental attempts to measure it, however, have yielded no result. Since 
Ohmic contact measurements cannot provide a result due to the constant elec- 
trochemical potential, 12,13 Bok and Klein 14 and later more precisely Morris and 
Brown 15 measured the Bernoulli potential at the surface via capacitive pickup. For 
the investigated Pb at 7K the thermodynamic corrections have been expected to 
be larger by a factor 30 than the first term in (11). Surprisingly the observed data 
agree only with the first term of (11). 

Why no signal of thermodynamic corrections is seen remained a puzzle for 
nearly 30 years. Recently we found the solution by a modification 16 of the Budd- 
Vannimenus theorem 17 

-*"»->- --»*!;(£)■ (12) 
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which states that the extension of the surface profile towards the surface, ip(0) has a 
step at the surface due to the surface dipoles which is given exclusively in terms of 
the electronic free energy with no regard to the potential inside. Writing the change 
of the free energy as Rickayzen, etp = — -§^Sf and Sf = n s \mv 2 , the potential at 
the surface becomes 

9 Sf 9„ n s 1 2 

e^surf = + n- = eip + — Sf = -mv . (13) 

On n on n n 2 

The surface dipoles compensate the thermodynamic corrections exactly. 16 

This is of course strictly valid only for homogeneous superconductors. Even the 
appearance of vortices in type-II superconductors requires a more detailed discus- 
sion. In order to describe such inhomogeneous situation we use the extension of the 
Landau-Ginzburg equation towards lower temperatures by Bardeen. 18 



2.3. Inhomogeneous superconductors within extended 
Ginzburg- Landau approach 

We consider the free energy as composed of the condensation energy the kinetic en- 
ergy, the electrostatic energy as well as the magnetic energy, !F[ip, A, n n ] = f drf = 
T s +-7 7 kin +J~c +^ 7 m- The condensation energy is constructed according to the two- 
fluid model of Gorter and Casimir, 19 where one describes two experimental facts, 
firstly that the fraction of superconducting electrons shows a temperature depen- 
dence of nearly w = 1 — ^ rts - " s w and secondly that the condensation energy is 
an expression containing the specific heat parameter 7 and the critical temperature 
as e con = \lT 2 ■ Eliminating the critical temperature one sees that the condensation 
energy must have the form e con = 4 jf_ ro in equilibrium. This relation should result 
from a variation of the free energy SF s jSvD = df s /dw = such that 

f s = U- e con w - ]pT 2 VT^. (14) 

The kinetic part is taken as a form proposed by Ginzburg and Landau, 20 who sug- 
gested to represent the superconducting density by a pseudo wave function \tp\ 2 = 7J f 
such that the kinetic energy can be given in Schrodinger form (m* = 2m e , e* = 2e) 

.Fkin= fdr^\(-ihV-e*A)i;\ 2 . (15) 

Variation of the total free energy with respect to ip* yields the Ginzburg-Landau 
equation (GL) for the wave function 

1 



2m 

with the potential 



-(-ifiV- e*A) 2 V + xV' = (16) 



9fs r,df s _ £con 1T 2 1 

X c#| 2 dn n n 2n L 2\w ' 

Y n 
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This extends the GL equation towards lower temperatures. Close to T c the potential 
takes the form \ ~~ * a + flW 2 an d the original GL equation is recovered. 

The GL equation was solved numerically 21 for Nb; the results are seen in figure 1. 
The vortices arrange themselves in the hexagonal structure of an Abrikosov-type 
lattice. The superconducting density is zero at the vortex centers and reaches nearly 
the nonmagnetic value at the borders of the vortex. The fact that n s is smaller there 
than its nonmagnetic value is a result of the nonlocal effects. The corresponding 
magnetic field reaches its maximum -B max at the vortex centers. In fact, it is higher 
than the applied field in the core which means that the superconductor compresses 
the magnetic field in vortices. 




2.4. Charge profile 

The numerical solution of the GL equation allows us to discuss the detailed charge 
density and potential profile towards the surface. For magnetic fields perpendicular 
to the surface we found a linearized solution of the GL equation (16) with respect 
to low magnetic field ip = y/n(l - (T/T c ) 4 )/2 + Sip with 



22 



where A = m* /fiQiie 2 is the London penetration depth indicating how far the mag- 
netic field penetrates the superconductor, £ = h 2 /2m*a is the GL coherence length 
in terms of a = 7T C (T — T c ) j n, and the ratio between both quantities is the GL pa- 
rameter k. One sees that the two length scales determine the profile corresponding 
to whether k is larger or smaller than \/2, i.e. type II or type I superconductors. The 
resulting charge profile results from the variation of the free energy with respect to 
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the normal density and reads 



1 2 ( n s t n n d\nT c \ 1 cv A « 



2 V n n <91nn / i + VI \ 1 - 

This result remains finite at the limiting case of k = \/2 and yields Rickayzen's 
result (11) for k — ► oo. 

With the help of the Poisson equation the linearized density profile reads 22 



(l — V n n <91nn y\ 

The surface value, p(x = 0), is always negative while the bulk value, p(x 3> 0), 
remains positive. This allows us to define the width of the surface by p(w) = with 
the result 

W= 9^TT- ( 21 ) 

For type-II superconductors, k > y/2, the surface charge is formed by the con- 
tribution on the scale of the GL coherence length £ while for n — > oo the width 
w — » In and the bulk charge extends on the scale of A. For type-I super- 
conductors, /t < \/2, the situation is reversed, the surface charge is formed by the 
contribution on the scale of the London penetration depth A and for n <C V2 the 
width w — > § In ^2 an d the bulk charge extends on the scale of £. 



3. Application of the Bernoulli potential 
3.1. Surface potential 

Summarizing, the electrostatic potential can leak out of a superconductor by three 
types of charges: (i) The bulk charge which describes the transfer of electrons from 
the inner to the outer regions of vortices creating a Coulomb force. This force has 
to balance the centrifugal force by the electrons rotating around the vortex center, 
the outward push of the magnetic field via the Lorentz force and the outward force 
coming from the fact that the energy of Cooper pairs is lower than the one of free 
electrons such that unpaired electrons in the vortex core are attracted towards the 
condensate around the core. 22 (ii) The surface dipole which cancels all contributions 
of pairing forces 23 resulting in an observable surface potential of 

e^o = - — • (22) 
n 

(iii) The surface charge 24 distributed on the scale of the Thomas-Fermi screening 
length Atf with ^,$o,A> Atf — > 0. 

Within the discussed Bardeen's extension of the GL equation the free electron 
energy has three components 
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which determines, according to the modified Budd-Vannimenus theorem (22), the 
surface potential. Near the critical temperature this becomes very simple e<f)o — 
^f3\4>\ 4 in terms of (3 = -fT 2 /2n 2 . Without surface dipolcs, the surface potential 
equals the internal potential 

where the first part are the inertial and Lorentz forces neglecting pairing 
forces. The approximation by Khomskii and Kusmartsev 25 adopted by Blatter, 26 
e0Bi = "^T^f W 2 1 takes only the pairing force into account. Before comparing the 
different approximations of the literature with our result, it is illustrative to discuss 
the magnetic field and temperature dependence. For thin layers and a B field close 
to B c2 the spatially averaged values of the superconducting fraction uj = , 2 can 

be given by general arguments 27 as (uj) = ^ and (uj 2 ) = with a number (3a 

characteristic for the vortex geometry and b = -g^-. Knowing our averaged surface 
potential we compare it now with Blatter's approximation 

( e ^o) = ^(l-t 2 ) 2 (l-6) 2 (25) 

where t = T/T c . Obviously the surface dipole strongly modifies the magnitude of 
the potential, in particular when the GL wave function has a small magnitude we 
have 4>q oc l^l 4 , while without dipoles 0bi oc \ip\ 2 . 

Possible experimental consequences for which the presented theory can be tested 
are at temperatures close the critical temperature, t — > 1, where l^l 2 oc 1 — t, 
therefore cf>o oc (1 — t) 2 while </>bi oc 1 — t and at magnetic fields close to the upper 
critical field, b — > 1, where l^l 2 oc 1 — 6 so that ^oc(l- b) 2 while B i oc 1 — b. 

The numerical comparison 24 for Nb is shown in figure 2. The internal potential 
and Blatter's result are similar. The full theory and the approximations of only iner- 
tial/Lorentz forces are much smaller. The surface dipole cancels a major part of the 
pairing forces. The full theory and the incrtial/Lorentz forces even result in different 
profiles and sign. This illustrates the delicate balance between the three contribu- 
tions in (23). It is not justified to neglect just one of them. Corresponding results 
in the literature should be critically reexamined, e.g. the vortex induced deforma- 
tion and the experimentally accessible effective mass of vortices show remarkably 
different results. 28 
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2. (a) Surface potential for different approximations at 1 



-7T- - 0.95, k = 0.78 and -g- = 
0.7818. The broken lines indicate the behavior along y directions different from the x directions 
due to the hexagonal structure, (b) An enlarged view of the approximation of only inertial and 
Lorentz forces compared to the present theory. 24 



3.2. Charged vortices probed by NMR 

Though the capacitive coupling measurements of the surface potential do not yield 
an access to thermodynamic corrections and the high precision measurements of 
the electrostatic potential leaking out of the surface remain to be performed, there 
is a direct attempt to measure the charge deeper in the bulk by NMR. Kumagai et 
al. 29 have measured the quadrupole resonance lines in the high-T c material YBCO. 
The polarization of the Cu atoms leads to a coupling of spins with the electrical 
field gradient and to a splitting of the quadrupole resonance 

z^ QR = E ±3/2 - E ±1/2 =A P + C (27) 

proportional to the charge density p. This allows to measure the charge of the 
vortices. If one advocates the picture that the positive charges are expelled from 
the vortex to the outer regions one would expect a BCS estimate of the measured 

charge per length Q « ^pf- f ^pj ( d'ln "/T ) • This mcans f° r the underdoped regime 
Q > and for the overdoped Q < 0, which is just the opposite to what is observed 
experimentally. Moreover the BCS estimate gives a charge of 10 -5 e per vortex 
while the experiment yielded a value three orders of magnitude larger. We solved 
this apparent contradiction with the help of the Bernoulli potential as follows. 

One has to consider that YBCO is built up from Cu chains and planes where 
in the latter the superconductivity and the vortices are located. The discussed 
Bernoulli potential therefore appears in the planes. On the other hand the charge 
neutrality requires that the exceeding charge has to be compensated in the Cu 
planes such that a charge transfer appears between planes and chains. Taking this 
and the two-dimensional geometry into account we have been able to explain the 
observed quadrupole line shift, 30 see figure 3. The experimental line width is com- 
parable to the spatial variation of the charge and magnetic field profile. Therefore 
the actual spatial profile of the magnetic field and charge density is important as 
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T/T, 

(a) (b) 

Fig. 3. (a) Spatial variation of the line shift V2 = E—1/2 — E_ 3 / 2 - The units are in vortex sizes, 
(b) Shift (Au) of r = 140kHz of a single crystal (dashed line), and after averaging over grain 
orientation (solid line) compared to the experimental values. 21 



presented in figure 3. 

4. Consequences to the kinetic theory 
4.1. Extended quasiparticle approximation 

We have seen that the concept of correlated density leads in a natural way to the 
Bernoulli potential, which is useful in explaining recent experiments and promising 
to access a material parameter in superconductors. The underlying noncquilibrium 
theory has to distinguish between parts of the Wigner distribution described by 
quasiparticle poles and the rest called off-shell parts. This off-shell motion can be 
eliminated from the kinetic equation, which requires to introduce an effective dis- 
tribution (the quasiparticle distribution /) from which the Wigner distribution fw 
can be constructed 

M/]-/*/!^ ((!-/>•£-/.*). < 28 » 

Here a > and cr < denote the selfenergies describing all correlations and e is the 
quasiparticle energy. This relation represents the extended quasiparticle picture 
derived for small scattering rates. 5,31-33 The limit of small scattering rates has been 
first introduced by Craig. 34 An inverse relation f[fw] has been constructed. 35 For 
equilibrium non-ideal plasmas this approximation has been employed by 36,37 and 
has been used under the name of the generalized Beth-Uhlenbeck approach by 38 
in nuclear matter for studies of the correlated density. The authors in 39 have used 
this approximation with the name "extended quasiparticle approximation" for the 
study of the mean removal energy and high-momenta tails of Wigncr's distribution. 
The non-equilibrium form has been derived finally as the modified Kadanoff and 
Baym ansatz. We will call this the extended quasiparticle approximation. 

This extended quasiparticle picture leads to balance equations which include 
explicit correlation parts analogous to the virial corrections. The firmly established 
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concept of the equilibrium virial expansion has been extended to noncquilibrium 
systems 5 although a number of attempts have been made to modify the Boltzmann 
equation so that its equilibrium limit would cover at least the second virial coeffi- 
cient. 41-43 The corrections to the Boltzmann equation have the form of gradients or 
nonlocal contributions to the scattering integral. The nonlocal quasiparticlc kinetic 
equation for the momentum-, space-, and time-dependent distribution function of 
particle a, fi = f a (k,r,t), derived within the non-equilibrium Green's function 
technique 5 ' 31 has the form of a Boltzmann equation with the quasiparticle energy 
si = s(k,r,t), 

^ll + ^l^Jj_ _ ^£1^/1 _ s y^ f d P d 1 

r)t flk r)r fir flk ^— ' 



(2ir) 5 h 7 



(l-/l)(l-/ 2 -)/ 3 "/4- - /l/ 2 ± (l-/ 3 ± )(l-/4 ± ) 



(29) 



and the spin-isospin etc. degeneracy s. The superscripts ± denote the signs of non- 
local corrections: f£ = f b (p,r±A 2 ,t), f% = f a (k - q± A K , r± A 3 , t ± A t ), and 
ff = fb(p+q±Ax, r±Ai, t±A t ). For the out-scattering part of (29) either the plus 
or minus signs can be chosen. 44 The scattering measure is given by the modulus 
of the scattering T-matrix V± = \T±\ 2 5(si + £2 — £3 — £4 ± 2A_e). All corrections 
A, describing the nonlocal and non-instant collision are given by derivatives of the 
scattering phase shift = Im \nT R (Q, k,p,q,t,r) 



A t = 



A 2 = 



S1+S2 



dp dq dk 



E1+E2 



A, 



ld<j> 
2~dt 



S1+S2 



dk 



(30) 



si+e 2 



A K = 



1 

2 dr 



dj)_ 
dk 



S1+S2 \ y / £l+£2 

The nonlocal kinetic equation (29) covers all quantum virial corrections on the bi- 
nary level and conserves density, momentum and energy including the corresponding 
two-particle correlated parts. 5 It requires no more computational power than solving 
the Boltzmann equation. 45,46 



4.2. Conservation Laws 

Neglecting all shifts, the time-invariant observables are the mean quasiparticle den- 
sity n qp , the mean momentum <2 qp , the mass current j qp , the mean energy £ qp and 
the stress tensor J°f in the form 

n qp - E /, 2 qp = E k /, j qp = E ff/> 

k k k /qi \ 

£ qp = + kmf)/ fe , J,f = E (kj f- + Sue) f %£ qp , 
k fe v y 
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where we abbreviated the mean-field part rj m f = / t R {u> + s p , k,p, 0)f p of the 
sclfcnergy. 

Taking now into account the binary correlations which we have reformulated 
into the shifts, we obtain from the non-local kinetic equation the modified balance 
equations 5 

g (n frcc + n corr) ^qp + jcorr) _ 

dt + d~r " u {6 > 

d(Qf + Qf rr ) , x d(Jff + Jl 



E fP 1 ' = (33) 



dt t-r* dr. 

d£ _ <9(£ qp + £ 



dt dt =°- (34) 

This shows that the conserving observables now consist of the sum of the quasipar- 
ticle parts (31) and the correlated parts which we can interpret as the parts coming 
from correlated pairs or molecules, 



n corr 


= JdVA t 


f°" = J dVA 3 


Qcorr 


= JdV^ 


A t £ co " - j dV ek ^ tp A t 






kjA 3i + Pj(A 4i - A 2 ») + qj(A 4i 



^3i)|, (35) 

where dV is the probability to form a pair per unit of time. This immediate inter- 
pretation of dV is obvious from the expressions (35). The density of pairs is given, 
if dP is multiplied by the lifetime of the molecule At. The energy, the momentum 
and mass current which are carried by the pairs result analogously. As one can 
see, the correlated part of the stress tensor (35) takes the form of a virial which is 
well known for the collision flux in dense gases. The corresponding momenta are 
connected to the corresponding offsets. The diagonal part of this stress tensor is 
of course the correlated part of the pressure. The proofs that these conservation 
laws indeed follow from the nonlocal kinetic equation and are consistent with the 
extended quasiparticle picture can be found in detail in the appendix G of. 5 



5. Conclusion 

The nonlocal kinetic equation unifying the achievements of transport theory in dense 
gases with the quantum transport of dense Fermi systems is presented. The quasi- 
particle drift of Landau's equation for the evolution of the quasiparticle distribution 
is connected with a dissipation governed by a nonlocal and non-instantaneous scat- 
tering integral in the spirit of Enskog corrections. 47 In this way quantum transport 
can be recast into a quasiclassical picture. The balance equations for the density, 
momentum and energy include quasiparticle contributions and the correlated two- 
particle contributions beyond the Landau theory. 48 

Compared with the Boltzmann equation, the presented form of virial corrections 
only slightly increases the numerical demands in implementations. Within INDRA 
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collaboration, we have performed numerical studies and compared them with the 
experimental data from GANIL. The nonlocal corrections with parameters found 
from collisions of two isolated nucleons 47 have been implemented into numerical 
simulations of heavy ion reactions in the non-relativistic regime. 45 ' Keeping all 
simulation parameters as in the local approximation, the temperature of mono- 
nucleon products of central reactions increases due to nonlocal effects towards ex- 
perimentally observed values, while the temperature of more complex clusters re- 
mains unchanged. 45 For the proton distribution in very peripheral reactions we have 
achieved an agreement between theoretical predictions and experimental data. 46 

As a consequence of the nonlocal kinetic theory, the correlated density appears 
which results in the Bernoulli potential in superconductors. This Bernoulli potential 
allows to access material parameters. Since at the surface large-scale cancellations 
appear, the charge of vortices as measured by NMR could be explained. We believe 
the correlated density to be a fruitful concept to describe strong correlations in 
interacting systems. 
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